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CN ' Abstract 

Q , Fractals are measurable metric sets with non-integer Hausdorff dimensions. If elec- 

\ trie and magnetic fields are defined on fractal and do not exist outside of fractal in 

Euclidean space, then we can use the fractional generalization of the integral Maxwell 
equations. The fractional integrals are considered as approximations of integrals on 
fractals. We prove that fractal can be described as a specific medium. 
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1 Introduction 

The theory of integrals and derivatives of non-integer order goes back to Leibniz, Liouville, 
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, Riemann, Grunwald, and Letnikov [H [2] . Fractional analysis has found many applications in 

recent studies in mechanics and physics. The interest in fractional integrals and derivatives 
has been growing continuously during the last few years because of numerous applications. 
In a short period of time the list of applications has been expanding. It includes chaotic 
O ■ dynamics [H H], physics of fractal and complex media O El El [8], [9], physical kinetics [HI 

flOl [TTl [T2] , plasma physics [IHl E] , astrophysics [I5] , long-range dissipation [161 [E] i non- 
^ \ Hamiltonian mechanics [iHl [19] , theory of long-range interactions [201 EH ES] . 

^ \ The new type of problem has increased rapidly in areas in which the fractal features of 

the process or the medium impose the necessity of using non-traditional tools. In order to 
use fractional derivatives and fractional integrals for media on fractal, we must use some 
continuous model [H]. We propose to describe the medium on fractal by a fractional con- 
tinuous model [H], where all characteristics and fields are defined everywhere in the volume 
but they follow some generalized equations, which are derived by using fractional integrals. 
In many problems the real fractal structure can be disregarded and the medium on fractal 
can be described by some fractional continuous mathematical model. The order of the frac- 
tional integral is equal to the fractal dimension. Fractional integrals can be considered as 
approximations of integrals on fractals [231 [21] • In Ref. [21], authors proved that integrals 
on a net of fractals can be approximated by fractional integrals. In Ref. [18], we proved 
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that fractional integrals can be considered as integrals over the space with non-integer di- 
mension up to numerical factor. This interpretation follows from the well-known formulas 
for dimensional regularization [29] . 

In Sec. 2, a brief review of Hausdorff measure, Hausdorff dimension and integration 
on fractals suggested to fix notation and provide a convenient reference. The connection 
integration on fractals and fractional integration is discussed. In Sec. 3, the fractional 
electrodynamics on fractals is considered. Fractional generalization of the integral Maxwell 
equations is suggested. Finally, a short conclusion is given in Sec. 4. 



2 Integration on fractal and fractional integration 
2.1 Hausdorff measure and Hausdorff dimension 

Fractals are measurable metric sets with fractal Hausdorff dimension. The main property 
of fractal is non-integer Hausdorff dimension. Let us consider a brief review of Hausdorff 
measure and Hausdorff dimension to fix notation and provide a convenient reference. 

Consider a measurable metric set {W, fin) with W C M". The elements of W are denoted 
by x,y,z,..., and represented by n-tuples of real numbers x = (xi, X2, such that W 
is embedded in M"^. The set W is restricted by the conditions: (1) W is closed; (2) W is 
unbounded; (3) W is regular (homogeneous, uniform) with its points randomly distributed. 

The diameter of a subset E C IV C M" is 

d{E) = diam{E) = sup{d{x,y) : x,y E E}, 

where d{x, y) is a metric function of two points: x and y G W. 

Let us consider a set {Ei} of subsets Ei such that diam{Ei) < e for all i E N, and 
W C Ui^i ^i- Then, we define 

aE^,D) = uj{D)[diam{E,)]'' = uj{D)[d{E,)]'' (1) 

for non-empty subsets Ei of W. The factor uj{D) depends on the geometry of Ei, used for 
covering W. If {Ei} is the set of all (closed or open) balls in W, then 

The Hausdorff dimension D of a subset E C W is defined [25l[26l[271[28] bv 

D = dimniE) = snp{d e R : finiE, d) = oo} = mf{d G M : finiE, d) = 0}. (3) 
The Hausdorff measure of a subset E cW is defined [23 ESI EZl EB] by 

oo 

fiH{E,D)=u{D) hm M\2[d{E,)]^. (4) 

d{E,)^0+ {E,} 

Note that fiH{XE,D) = X^fiH{E,D), where A > 0, and = {Ax, x G E}. 
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2.2 Function and integrals on fractal 

Let us consider the functions on W: 

oo 

f{x) = '^(^tXEA^)^ (5) 

i=l 

where xe is the characteristic function of E: 

XEix) 

The Lebesgue-Stieltjes integral for ([5]) is defined by 

„ oo 

/ fd^I = y2^^^ME^). (6) 

Jw 



1 if X e E, 
if X ^ E 



Therefore 



/ f{x)dfiH{x)= hm y2f{x,)aE^,D)=uj{D) hm V /(x^) 



(7) 

Er 

It is possible to divide into parallelepipeds 

Ei^...i„ = {{xi, ...,Xn) eW : Xj = [ij - l)Axj + a^, < < Axj, j = 1, ...,n}. (8) 
Then 

n n 

The range of integration W can be parametrized by polar coordinates with r = d{x, 0) and 
angle Q. Then ii^^.n can be thought of as spherically symmetric covering around a center at 
the origin. In the limit, the function C,{Er^n,D) gives 

diiH{r,n)= lim i{Ern,D)=dVl^'^r^~^dr. (10) 

Let us consider f{x) that is symmetric with respect to some center xq G W , i.e. f{x) = 
const for all x such that Xq) = r for arbitrary values of r. Then the transformation 

W^W: x^x' = x-Xo (11) 
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can be performed to shift the center of symmetry. Since W is not a hnear space, (JTTl) need 
not be a map of W onto itself. The map (JTTl) is measure-preserving. Then the integral over 
a D-dimensional metric space is 

This integral is known in the theory of the fractional calculus [1]. The right Riemann- 
Liouville fractional integral is 

/-/(^) = ^^y^ {x-z)^-'f{x)dx. (13) 

Equation ([T2|) is reproduced by 



r 27r^/^r(D) 



/(O). (14) 

Relation ( |T4|) connects the integral on fractal with integral of fractional order. This result 
permits to apply different tools of the fractional calculus [l] for the fractal medium. As a 
result, the fractional integral can be considered as an integral on fractal (fractional Hausdorff 
dimension set) up to the numerical factor T{D / 2) / [271^ ^"^T^D)]. 

Note that the interpretation of fractional integration is connected with fractional dimen- 
sion [18] . This interpretation follows from the well-known formulas for dimensional regular- 
ization [291 [30] . The fractional integral can be considered as an integral in the fractional 
dimension space up 

to the numerical factor r(i:)/2)/[27r^/2rp)]_ Ref. [23] was proved 
that the fractal space-time approach is technically identical with dimensional regularization. 
The integral defined in ([7j) satisfies the translational invariance property: 

f{x + Xo)dfiH{x) = / f{x)dfiHix) (15) 
w Jw 

since dfi^ix — Xq) = dfi^ix) as a consequence of homogeneity. The integral ([7]) satisfies the 
scaling property: 

' f{Xx)dfiH{x) = A-^ / f{x)dfiH{x) (16) 
w Jw 



since diJ,H{x/X) = A djj,H{x). 



2.3 Multi-variable integration on fractal 

The integral in (IT^ is defined for a single variable but not multiple variables. It is only 
useful for integrating spherically symmetric functions. We consider multiple variables by 
using the product spaces and product measures. 
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Let us consider a collection of n = 3 measure spaces {Wk, fJ^k, D) with k = 1,2,3, and 
form a Cartesian product of the sets Wk producing the space W = Wi x W2 x W3. The 
definition of product measures and application of Fubini's theorem provides a measure for 
the product set W = Wi x W2 y< W3 as 

(/il X /i2 X fis){W) = fIl{W,)^l2iW2)l^3{W,). (17) 

The integration over a function / on the product space is 

/(r)(i(/ii X /i2 X /is) = / / f{xi,X2,X3)dfii{xi)dH2{x2)dl^3{x3). (18) 



In this form, the single- variable measure from f[T^ may be used for each coordinate Xk, which 
has an associated dimension a^: 

27i-0!fc/2 

diik{xk) = —, — —-\xk\"''~'^dxk, k = 1,2,3. (19) 
r(afc/2) 

The total dimension of W = Wi x W2 x W3 is D = ai + a2 + a^. 

Let us reproduce the result for the single- variable integration (1121) . from the product 
space Wi X W2 x W3. We take a spherically symmetric function /(r) = f{xi,X2,xs) = f{r), 
where = (xi)^ + (^2)^ + (^3)^, and to perform the integration in spherical coordinates 
(r, 0, 9), we use 

rs,n^--xcos--.dx=^'f/^'^<f>, (20) 
io 2r((^ + ,.)/2)' 

where > 0, z/ > 0. Then Eq. f|T8|) becomes 

dni{xi)diJ,2ix2)diJ,3{x3)f{r) = y(^j:)/2) J ^^^^^ 

This equation describes integration over a spherically symmetric function in a D-dimensional 
space and reproduces the result ffT^ . 



3 Electrodynamics on fractal 
3.1 Electric charge of fractal set 

Let us consider the electric charge that is distributed on the measurable metric set W with 
the fractional Hausdorff dimension D. Suppose that the density of charge distribution is 
described by the function p{r,t). In this case, the total charge is defined by 

Qd{W) = p{r,t)dVD, dVo = dpi{xi)dp2{x2)dp3{x3) = C3{D,Y)dV3, (22) 
Jw 
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where dVs = dxdydz for Cartesian coordinates, dimniW) = D = ai + a2 + as, and 

87r^/2|3;|^l-l|^|^2-l|^|a3-l 

C3p, r) = . (23) 

r(ai)r(a2)r(a3) 

As a result, we get Riemann-Liouville fractional integral |l] up to numerical factor Svr^/^. 
Note that the final equations that relate the physical variables have the form that are in- 
dependent of numerical factor in the function C3{D,t). However, the dependence of r is 
important to these equations. 

Equation (122|) describes the charge that is distributed in the volume and has the frac- 
tal dimension D by fractional integrals. There are many different definitions of fractional 
integrals [1]. For the Riemann-Liouville fractional integral, the function C3{D,r) is 

cs{D, r) = ' ' ; '/ , , (24) 

where x, y, z are Cartesian's coordinates, and D = ai + a2 + as, < D < 3. Note that 
for D = 2, we have the distribution in the volume. In general, this case is not equivalent to 
the distribution on the two-dimensional surface. For p(r) = p(|r|), we can use the fractional 
integrals with 

C3(Ar)^^9|fw- (25, 

If we consider the ball region W = {r : |r| < R}, and stationary spherically symmetric 
distribution of charged particles (p(r, t) = p(r)), then 



For the homogeneous case, p{r,t) = po, and 

QoiR) - 4.p„29|p^ ^ R^. 

The distribution of charged particles is homogeneous if all regions W and W with equal 
volumes VoiW) = Vd(W') have the same total charges on these regions Qd{W) = QniW'). 

For charged particles that are distributed with a constant density over a fractal with 
Hausdorff dimension D, the electric charge Q satisfies the scaling law Q{R) ~ R^, whereas 
for a regular n-dimensional Euclidean object we have Q{R) ~ M". 

3.2 Electric current for fractal 

For charged particles with density p{r,t) flowing with velocity u = u(r, t), the resulting 
current density J(r, t) is 

J(r,t) =p(r,t)u(r,t). 
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The electric current I{S) is defined as the flux of electric charge. Measuring the fleld J(r, t) 
passing through a surface 5* = dW gives 

I{S) = ^jiS)= [ (J(r,t),f/S2), (26) 

where dS^ = dS2ii is a differential unit of area pointing perpendicular to the surface S, and 
the vector n = n^efc is a vector of normal. The fractional generalization of (126!) is 



I{S)= I (J(r,t),rfS,), (27) 

where 



dS, = C2{d,r)dS2, C2{d,r) = :^^\r\'^-\ (28) 

Note that C2(2,r) = 1 for d = 2. The boundary dW has the dimension d. In general, the 
dimension d is not equal to 2 and is not equal to {D — 1). 

3.3 Charge conservation for fractal 

The electric charge has a fundamental property established by numerous experiments: the 
velocity of charge change in region W bounded by the surface S = dW is equal to the flux 
of charge through this surface. This is known as the law of charge conservation: 

or, in the form 



d 
di 

In particular, when the surface S = dW is flxed, we can write 



/ p{r,t)dVD = - (f (J(r,t),rfS,). (29) 
Jw Jaw 



Using the fractional generalization of the Gauss's theorem (see the Appendix), we get 

/ (J(r,t),rfSrf)= / c-^\D,v)^U{d,v)UT,t))dVD. (31) 
Jaw Jw oxk V / 

Substitution of Eqs. ^ and (Ell) into Eq. ([291) gives 



C3(Ar)^^ + ^(c2(rf,r)J,(r,t)) =0. (32) 

As a result, we obtain the law of charge conservation in differential form (l32l) . This equation 
can be considered as a continuity equation for fractal. 
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3.4 Electric field and Coulomb's law 

For a continuous stationary distribution p(r'), the electric field at a point r is defined by 

where Eq is a fundamental constant called the permittivity of free space. For Cartesian's 
coordinates dV^ = dx'dy'dz' . The fractional generalization of ( l33l) is 

where dV[) = c^^D, r')dV^. Equation ( l34l) can be considered as Coulomb's law for a stationary 
distribution of electric charges on fractal. 

3.5 Gauss's law for fractal 

The Gauss's law tells us that the total flux ^e{S) of the electric field E through a closed 
surface S = dW is proportional to the total electric charge Q{W) inside the surface: 

^E{dW) = -Q{W). (35) 

The electric flux for the surface S = dW is 



^EiS)= [ (E,rfS2) 

Js 



where E(r, t) is the electric field vector, and dS2 is a differential unit of area pointing per- 
pendicular to the surface S. 

For the distribution on fractal, the Gauss's law fl35l) states 

E,dS2) = - f p{r,t)dVD, (36) 
S ^0 Jw 

where p(r, t) is the density of electric charge that is distributed on fractal, dVo = c^i^D, r)dV3, 
and Eq is the permittivity of free space. 

If p{r,t) = p{r), and W = {r : |r| < R}, then 

Q{W)=4nj^ p(r)cs{D,rydr = 4n^,^^j^ piry^^'^dr. (37) 

For the sphere S = dW = {r : |r| = R}, 

(^E{dW) =A'kR^E{R). (38) 
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Substituting (EZD and (JMl) in (I35l), we get 



-SSI 

For homogeneous (p(r) = p) distribution, 

^ ' ^eoDV{D/2) 
3.6 Magnetic field and Biot-Savart law 

The Biot-Savart law relates magnetic fields to the currents, which are their sources. For a 
continuous distribution, the law is 

BW^£/„,^^;^-3, (3.) 

where [ , ] is a vector product, J is the current density, /io is the permeability of free space. 
The fractional generalization of Eq. (|39ll is 



Bfr) = ^ 



^ / ^^^^^dV^,. (40) 



This equation is Biot-Savart law written for a steady current with fractal distribution of 
electric charges. The law ( HOj) can be used to find the magnetic field produced by any 
distribution of steady currents on fractal. 



3.7 Ampere's law for fractal 

The magnetic field in space around an electric current is proportional to the electric current, 
which serves as its source. In the case of static electric field, the line integral of the magnetic 
field around a closed loop is proportional to the electric current fiowing through the loop. 
The Ampere's law is equivalent to the steady state of the integral Maxwell equation in free 
space, and relates the spatially varying magnetic field B(r) to the current density J(r). 

Note that, as mentioned in Ref. |l3j, Liouville, who was one of the pioneers in develop- 
ing fractional calculus, was inspired by the problem of fundamental force law in Ampere's 
electrodynamics and used fractional differential equation in that problem. 

The Ampere's law states that the line integral of the magnetic field B along the closed 
path L around a current given in MKS by 

jf(B,rfl)=/io/(5), 
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where dl is the differential length element. For the distribution of particles on the fractal, 
I{S) is defined by (1271) . For the cylindrically symmetric distribution, 



I{S) = 2n r J{r)c2{d,v)rdr = Att-^^ r J{r)r''-^dT, 
Jo r(rf/2) Jo 

where we use C2{d,r) from Eq. (l28l) . For the circle L = dW = {r : |r| = R}, we get 

(B,(il) = 27rR B{R). 



As a result. 

For the homogeneous case, J(r) = Jq and 



3.8 Fractional integral Maxwell equations 

Let us consider the fractional integral Maxwell equations. The Maxwell equations are the set 
of fundamental equations for electric and magnetic fields. The equations can be expressed 
in integral form known as Gauss's law, Faraday's law, the absence of magnetic monopoles, 
and Ampere's law with displacement current. In MKS, these become 

/(E,rfS2) = - / pdVo, 
Js ^0 Jw 

j((E,rfIi) = -^^(B,dS2), 
jf(B,rfS2)=0, 

y"(B,dli) =/io j (J, c/Sd) + eo/io^ j (E,c/S2). 

Let us consider the fields that are defined on fractal only. The hydrodynamic and thermo- 
dynamics fields can be defined in the fractal media p[|9]. Suppose that the electromagnetic 
field can be defined on fractal as an approximation of some real case with fractal medium. 
If the electric field E(r,t) and magnetic fields B(r, t) can be defined on fractal and do not 
exist outside of fractal in Euclidean space M^, then we must use the fractional generalization 
of the integral Maxwell equations in the form: 

(f{E,dS,) = - [ pdVn, (41) 
Js ^0 Jw 
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;E,rfg = -^ /(B,ds,), (42) 



;B,dS,) = 0, (43) 



^ (B, d\^) = /io ^(J, dSd) + sofio-^ J (E, dS, 



(44) 



The fractional integrals are considered as approximation of integrals on fractals [231 [2] 

Using the fractional generalization of Stokes's and Gauss's theorems (see the Appendix), 
we can rewrite Eqs. f HTj) - (jH]) in the form 



/ c^\D,r)div{c2id,r)E)dVD = — [ pdVo, (45) 

j^c,\d,r){curl{c,{^,r)E),dS,) = -^^(B,rfSrf), (46) 

/ c^\D, r)div{c2id, r)B)dVd = 0, 
Jw 

J {d, r) {curl {ci{-f, r)B), dS d) = jJ^o J {J,dSd) + SojJo^ J (E,dSd). 
As a result, we obtain 

div(c2{d,r)E) = -C3(D,r)p, (49) 



(47) 
(48) 



CMr/(ci(7,r)E) = -C2(rf,r)-B, (50) 

div (^C2id, r)B^ =0, (51) 

/ \ (9E 

CMr/^Ci(7,r)Bj = /ioC2(rf, r) J + £:oAioC2('i, r)— . (52) 

Note that the law of absence of magnetic monopoles for the fractal leads us to the equation 
div{c2{d,r)'B) = 0. It can be rewritten as 

divB = —(B, gradc2{d,r)). 

In the general case {d ^ 2), the vector grad {c2{d,r)) is not equal to zero and the magnetic 
field satisfies divB 7^ 0. If d = 2, we have dw{B) ^ only for non-solenoidal field B. 
Therefore the magnetic field on the fractal is similar to the non-solenoidal field. As a result, 
the field on fractal can be considered as a field with some "fractional magnetic monopole" 
Qm ~ (B, Vcs). 
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3.9 Fractal as effective medium 

The Maxwell equations P9l) - fl52l) on fractal can be considered as the equations for medium 



curl I H 



eff 
Pfreey 


(53) 


dt ' 


(54) 


1=0, 


(55) 


dt ■ 


(56) 



The effective Maxwell equations (!53l) -( j56l) prove that fractal creates some polarization 
and magnetization. In the equations, we use some effective fields 

E'f^ir, t) = ci(7, r)E(r, t), B^^^(r, t) = c^^d, r)B(r, t). (57) 

The fields W^f and B^^^ mean that electromagnetic fields E and B of particles are changed 
by fractal. 

Equations (I53l) - fl56l) have the effective free charge and current densities: 

ptLi^^t) = Cs{D,r)p{r,t), r"{r,t) = C2(rf, r) J(r, t). (58) 

We can interpret the existence of Pj'|!fg and J'^^^ in the equations as an effect of change of 
the free electric charge and current densities by fractal. This change exists in addition to 
the effect of appearance the dipole charges and polarization or magnetization currents. The 
fractal can be considered as a medium that has the electrical and magnetic permittivities in 
the form 

e = C2((i,r)c7^(7,r), = C2(rf, r)q ^(7, r). (59) 
The fields D and H are related to E"^-^-^ and B^-^-^ by the usual equations: 

D = eeoE^f^, H = — B"-^^. (60) 

Note that the continuity equation (152]) for fractal can be presented by 

dp^^f{r,t) 



dt 



+ dw(^3'ff{r,t)^ =0. (61) 



As a result, the fractal can be considered as a specific medium that changes the fields, 
free charges and currents in addition to the creation of polarization and magnetization. 
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4 Conclusion 

Fractals are measurable metric sets with non-integer Hausdorff dimensions. We consider 
the electric and magnetic fields that are defined on fractal and does not exist outside of 
fractal in Euclidean space. For charged particles that are distributed with a constant density 
over a fractal with Hausdorff dimension D, the electric charge Q satisfies the scaling law 
Q{R) ~ 1 whereas for a regular n-dimensional Euclidean object we have Q{R) ~ M"^. 
This property can be used to measure the fractal Hausdorff dimension D. 

The fractional integrals can be used to describe electromagnetic fields on fractals. These 
integrals are considered as approximations of integrals on fractals. The fractional generaliza- 
tions of integral Maxwell equations for fractal set are derived. The magnetic field on fractal 
can be considered as a field with some "fractional magnetic monopole". 

We can interpret the equations for electromagnetic fields on fractal as an effect of creation 
of some polarization and magnetization by fractal. Moreover, the electromagnetic fields are 
also changed by fractal. From the generalized Maxwell equations, we can see the effect of 
change of the free electric charge and current densities by fractal. This change exists in 
addition to the effect of appearance the dipole charges and polarization or magnetization 
currents. The electrical permittivity e and the magnetic permittivity yU of fractal are defined 
by the Hausdorff measure and dimension of fractal. 



Appendix: Fractional Gauss's theorem 



Let us derive the fractional generalization of the Gauss's theorem 




(62) 



where the vector J(r,t) = J^e^ is a field, and 



div{J) = dJ/dr = dJk/dxk- 



Here, we mean the sum on the repeated index k from 1 to 3 



Using 



dSd = C2{d, r)dS: 



•2, 



C2(rf,r) 



r(rf/2) 



r 



we get 




Note that C2(2, r) = 1 for c/ = 2. Using ([62]), we get 
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The relation dVs — c^^{D,r)dVD allows us to derive the fractional generalization of the 
Gauss's theorem: 

/ {J{r,t),dSa)= [ c^\D,r)dtv(c2{d,r)J{r,t)) dVo. 
Jaw Jw ^ ' 

Analogously, we can get the fractional generalization of Stokes's theorem in the form 

/(E,rfg= / C2-^(ci,r)(cMr/(ci(7,r)E),(iSd), 
Jl Js 

where 

2'-T(l/2) 
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